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ABSTRACT 



o 

■ We have obtained the average 1-500 keV spectrum of radio-quiet Seyfert Is using data 

\Q ; from EXOSAT, Ginga, HEAO-1, and GRO OSSE. The spectral fit to the combined 

Q\ ■ average EXOSAT and OSSE data is fully consistent with that for Ginga and OSSE, 

""^5 ' confirming results from an earlier Ginga/OSSF, sample. The average spectrum is well- 

fitted by a power-law X-ray continuum with an energy spectral index of a ~ 0.9 mod- 
erately absorbed by an ionized medium and with a Compton reflection component. A 
O ' high-energy cutoff (or a break) in the the power-law component at a few hundred keV 

or more is required by the data. We also show that the corresponding average spectrum 
from HEAO-1 Al and A4 is fully compatible with that obtained from EXOSAT, Ginga 
and OSSE. These results confirm that the apparent discrepancy between the results 
of Ginga (with a ~ 0.9) and the previous results of EXOSAT and HEAO-1 (with 
a ~ 0.7) is indeed due to ionized absorption and Compton reflection first taken into 
r> | account for Ginga but not for the previous missions. Also, our results confirm that the 

Seyfert- 1 spectra are on average cut off in 7-rays at energies of at least a few hundred 
keV, not at ~ 40 keV (as suggested earlier by OSSE data alone). The average spectrum 
is compatible with emission from either an optically-thin relativistic thermal plasma in 
a disk corona, or with a nonthermal plasma with a power-law injection of relativistic 
electrons. 

Key words: galaxies: active - galaxies: Seyfert - X-rays: galaxies - gamma-rays: 
observations - gamma-rays: theory - accretion, accretion disks 



1 INTRODUCTION 

The average X-ray/7-ray (hereafter X7) spectra of Seyfert 
Is and 2s observed by both Ginga and GRO OSSE have 
recently been obtained by Zdziarski et al. (1995, hereafter 
Z95). The main result of that study for Seyfert Is is that 
their average spectral high-energy cutoff is around several 
hundred keV. This is thus similar to the cutoff of IC 4329A, 
a bright Seyfert 1 (Madejski et al. 1995; Zdziarski et al. 1994, 
hereafter Z94). Z95 also found that the average spectrum of 
radio-quiet (hereafter RQ) Seyfert Is has the energy spectral 



index of a ~ 0.9 as well as a Compton-reflection component, 
as previously obtained for Ginga spectra alone (Pounds et 
al. 1990; Nandra & Pounds 1994, hereafter NP94). However, 
since the 23 Ginga and 8 OSSE observations of RQ Seyfert Is 
used in Z95 are not simultaneous and that sample consists of 
only 4 objects, there is a clear need to confirm those results 
using enlarged and independent samples. 

Here, we test and confirm the results of Z95 using spec- 
tra of RQ Seyfert Is from EXOSAT, HEAO-1, Ginga, and 
OSSE. [Discussion of the average X7 properties of radio- 
loud Seyferts is given in Wozniak et al. (1996).] We obtain 
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the average spectrum of RQ Seyfert Is detected by both 
EXOSAT and OSSE, a sample which consists of 7 objects. 
The objects were observed 41 times by EXOSAT and 18 
times by OSSE. Averaging this large numbers of observa- 
tions is expected to compensate for the lack of simultaneity 
of the observations. We also obtain the corresponding aver- 
age spectrum from HEAO-1 Al and A4. Furthermore, we 
analyze the Ginga/ OSSE sample enlarged by new OSSE ob- 
servations of NGC 5548 and with addition of NGC 7469. 

After presenting the spectra, we study physical pro- 
cesses that can be responsible for the observed X7 emission. 
We consider both thermal and nonthermal plasmas. 



2 THE DATA 

We use EXOSAT spectra from the HEASARC archive 
with the quality flag 3 or higher (which indicates obser- 
vations with relatively reliable background subtraction) of 
RQ Seyfert Is detected by OSSE. We exclude the Seyfert 
Is brightest in X-rays, NGC 4151 and IC 4329A, in order 
to avoid their dominance of the co-added spectrum. The us- 
able EXOSAT energy range is from 1.2 keV to 8 keV (chan- 
nels 6-31). The spectra above 8 keV suffer from relatively 
inaccurate global background subtraction [as opposed to lo- 
cal background subtraction used by, e.g., Turner & Pounds 
(1989)]. The individual spectra include a 1 per cent system- 
atic error. The spectra from EXOSAT (as well as Ginga 
and HEAO-1) are co-added with the weights corresponding 
to the length of time of each observation. Both the counts 
and the response matrices for each EXOSAT observation 
are added using a procedure (specially written for EXOSAT 
data) in the FTOOLS data processing package. 

The OSSE spectra take into account an estimated sys- 
tematic error correction to the spectra [see, e.g., Zdziarski, 
Johnson & Magdziarz (1996) for details]. We use the OSSE 
response matrix as revised in 1995, which results in the 50- 
60 keV fluxes about 20 per cent higher than in the earlier 
response (used, e.g., in Z95). 

The EXOSAT/OSSE sample consists of 7 RQ Seyfert 
Is [41, 18] (the numbers in brackets give the number of EX- 
OSAT and OSSE observations, respectively): MCG -6-30-15 
[4, 2], Mrk 509 [3, 2], NGC 3783 [4, 1], NGC 5548 [11, 6], 
MCG 8-11-11 [10, 2], ESO 141-55 [2, 2], and NGC 7469 [7, 
3]. 

We also use the average spectrum in the 2-180 keV 
range from observations by HEAO-1 Al and A4 (Wood et 
al. 1984; Rothschild et al. 1983). [We do not use HEAO- 
1 A2 data because the normalization of archival spectra is 
not known (Weaver, Arnaud & Mushotzky 1995).] The A4 
spectra of individual objects have been recreated using the 
current version of the instrument software. In the 2-10 keV 
energy range, we use the HEAO-1 Al data as published by 
Wood et al. (1984). Since that detector provides no spectral 
information and only instrumental counts are given in Wood 
et al. (1984), we have obtained the 2-10 keV fluxes using the 
counts-to-flux conversion as given for the Crab. In order to 
get an estimate of the 2-10 keV average spectrum, we use 
the average overall slope in that range of a — 0.7 (see Section 
[I] below). 

Our Ginga data are the same as those for Seyfert Is 
in Z95 except for the addition of NGC 7469 (NP94). The 



Ginga/ OSSE sample thus consists of 5 RQ Seyfert Is [25, 
14] (the numbers in brackets give the number of Ginga and 
OSSE observations, respectively): MCG -6-30-15 [4, 2], Mrk 
509 [4, 2], NGC 3783 [1, 1], NGC 5548 [14, 6], and NGC 7469 
[2, 3]. We have obtained the average spectra for both the 
top-layer and the mid-layer Ginga data (Turner et al. 1989). 
However, we have found that the mid-layer spectrum above 
10 keV is systematically softer than the corresponding top- 
layer spectrum. Since the mid-layer calibration is much more 
uncertain than that of the top layer (Turner et al. 1989), 
we use in this paper only the top-layer data. We use the 
energy range of 1.7-18 keV (channels 4-31), for which the 
instrumental background subtraction is accurate. As in Z95, 
a 0.5 per cent systematic error correction has been applied 
to the co-added Ginga spectrum. 



3 MODEL 

In our fits, we use the xspec spectral fitting package ver- 
sion 9.0 (Arnaud 1996). As in Z95, we model the underlying 
continuum as a power law with an exponential cutoff at an 
energy, E c . Z95 found that ionized low-energy absorption is 
necessary to explain the average Ginga spectrum of Seyfert 
Is. Thus, we model absorption as due to an ionized medium 
with the abundances from Anders & Ebihara (1982) and the 
ion opacities of Reilman & Manson (1979) except for the Fe 
K-edge energies, for which results of Kaastra & Mewe (1993) 
were used. The absorber temperature is assumed to equal 
T = 10 5 K (Krolik & Kallman 1984), its column density is 
./Vh, and an ionization parameter is defined by £ = L/(nr 2 ). 
Here L is the 5 eV-20 keV luminosity in a power law spec- 
trum with the average energy index of 0.7, and n is the 
density of the absorber located at distance r from the illu- 
minating source. Model parameters are given at the average 
redshift for each sample, (z). In addition to the ionized ab- 
sorber, we include a fixed neutral absorber at z — with 
Ah,g equal to the average Galactic value for each sample. 

The underlying continuum irradiates cold matter, e.g., 
an accretion disk, and gives rise to a Compton-reflection 
spectral component (Lightman & White 1988). Differently 
from Z95, who used the reflection spectrum averaged over 
all the angles of the reflected photons (Lightman & White 
1988), we assume a viewing angle of 30° [corresponding to 
an orientation close to face-on expected in type-1 AGNs, 
Antonucci (1993)]. We use the angle-dependent reflection 
Green's functions of Magdziarz & Zdziarski (1995). The lu- 
minosity intercepted by the reflecting medium equals R (the 
relative contribution of reflection) times the luminosity emit- 
ted outward, i.e., 2irR gives the solid angle subtended by the 
absorber as seen from the source of the incident radiation. 
The continuum reflection is accompanied by emission of a 
fluorescent Fe Ka line, which we model here as a Gaussian 
at an energy _Ef c and the width CFe- 



4 RESULTS 

Wc first study the Ginga/ OSSE average spectrum, using the 
sample enlarged with respect to that in Z95 (see Section ^). 
The data and the best-fit model are shown in Fig. la. The fit 
parameters are given in Table 1. We obtain a ~ 0.90 ± 0.05 
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Figure 1. The average spectrum (crosses) of RQ Seyfert 1 galax- 
ies from Gin 5 a/OSSE (a) and EXOSAT/OSSE (b). The upper 
limits here and in figures below are 2-<r. The dashed curves rep- 
resent power-law spectra with exponential cutoffs and bound-free 
absorbed at low energies, and the dotted curves represent Comp- 
ton reflection including the Fe K« line. The solid curves give the 
sum. 



and the e- folding energy between ~ 0.4 MeV and 2.7 MeV, 
which is in agreement with Z95. The relative contribution 
of reflection, 7? ~ 0.76 ± 0.15 is about 2/3 of that found by 
Z95, which is explained by the angle-dependent reflection 
spectrum at 30° both having a higher normalization with 
respect to the underlying continuum and being harder in 
the ~ 10-30 keV range than the angle-averaged reflection 
spectrum (Magdziarz & Zdziarski 1995). [Note that this ef- 
fect will also reduce large values of angle-averaged R ob- 
tained by Weaver et al. (1995) for Seyfert Is observed by 
HEAO-1 A2.] Our values of a and R are in agreement with 
the average values for Ginga of NP94, a ~ 0.95 with a dis- 
persion of 0.15 and R ~ 0.5-0.7. The equivalent width of 
the Kq line is 1201^ which agrees with the value ex- 
pected from fluorescence in the reflecting medium at our R 
(George & Fabian 1991). We have also confirmed that our 
results are only weakly dependent on the (relatively uncer- 
tain) elemental abundances. When the abundances of An- 
ders & Ebihara (1982) are replaced by those of Anders & 
Grevesse (1989) (with higher abundances of metals, in par- 
ticular with about 40 per cent more Fe), our results change 
only slightly. E.g., marginally more Compton reflection is 
obtained, R = 0.8012:1! (x 2 = 71.8/72 dof). 

Then we consider the EXOSAT/OSSE sample, see Fig. 



16. Since the usable EX OS AT spectrum extends only to 8 
keV, Compton reflection is not constrained. Therefore, we 
fix R at the value obtained from the Ginga/OSSE average. 
The equivalent width of the Kq line obtained, IOO^jq eV, 
is consistent with this assumption. (The best-fit line energy 
is 5 per cent less than that for Ginga, which is due to a gain 
inaccuracy of EXOSAT.) Furthermore, we find that there is 
an apparent soft X-ray excess present in the average spec- 
trum below 2 keV, a feature common in Seyfert 1 spectra 
(Wilkes & Elvis 1987; Turner & Pounds 1988). Rather than 
add a separate spectral component (which would be poorly 
constrained by our data), we fit only the spectrum above 2 
keV. We show, however, the spectrum below 2 keV in Fig. 
16. 

We find that the fit parameters of the EXOSAT/OSSE 
average spectrum are virtually identical to those for the 
Ginga/OSSE average. The two spectra are compared in Fig. 
2. We thus confirm that the average value of the e-folding 
energy in the spectra of Seyfert Is is several hundred keV 
rather than ~ 40 keV (Johnson et al. 1994). The latter value 
was obtained using narrow-band OSSE data only, neglecting 
Compton reflection, using a spectral model implying X-ray 
power laws much harder than observed, as well as with the 
old OSSE response. Those factors fully explain the discrep- 
ancy. 

We point out that the OSSE spectra of individual 
Seyferts appear more uniform (Johnson et al. 1994) than 
their X-ray spectra, for which the 1-a dispersion of a is 

0. 15 ± 0.04 (NP94). If indeed AGNs with different a have 
similar 7-ray spectra, there will be a positive correlation of 
E c with a, rather than a constancy of E c among objects 
with different a (because harder X-ray spectra will need to 
be cut off stronger than softer ones in order to yield sim- 
ilar spectra above 50 keV). The range of E c among indi- 
vidual AGNs will be then larger than that given in Table 

1, which range corresponds to our average X-ray spectrum, 
with a = 0.90 ± 0.05. This appears to be confirmed for 
NGC 4151, a bright Seyfert with a hard X-ray spectrum. 
Zdziarski et al. (1996) have found no statistical difference 
between the shape of the spectra of NGC 4151 from four 
OSSE observations during 1991-93 and the average OSSE 
spectrum of the Seyfert Is observed by EXOSAT. However, 
the e-folding energy for NGC 4151 is E c ~ 150 keV (below 
the range of E c in Table 1) as a consequence of a ~ 0.7 [for 
the cut-off power-law model with reflection fitted to 1991 
June Ginga/OSSE observation (Zdziarski et al. 1996)]. 

The agreement between the average EXOSAT and 
Ginga spectra confirms that the average value of a ~ 0.7 
for Seyferts found based on EXOSAT observations (Turner 
& Pounds 1989) is indeed an artifact of assuming neutral ab- 
sorption and not including Compton reflection. When both 
effects are included, the average a ~ 0.9 for both EXOSAT 
and Ginga. 

A similar average spectral index of a ~ 0.6-0.7 was also 
obtained by HEAO-1 (Rothschild et al. 1983; Mushotzky 
1984). Here, we obtained the average HEAO-1 A1/A4 spec- 
trum for the objects in the EXOSAT/OSSE sample (see Sec- 
tion |i|), shown as dotted symbols in Fig. 2. We see that apart 
from a small difference in the normalization the two spec- 
tra are fully consistent with each other. Thus, the spectra of 
Seyfert Is from EXOSAT and HEAO-1 are compatible with 
an intrinsic power law with a ~ 0.9. Our average HEAO-1 
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Table 1. Parameters of the spectral fits to the average RQ Scyfcrt 1 spectrum. The Ginga/ OSSE and EXOS AT/ OSSE samples are 
denoted by 'GO' and 'EO', respectively. The symbols arc explained in the text. Parameters with no error ranges are fixed. are in 
units of 10 21 cm~ 2 , E c , E Fc , and <xp e are in keV. Errors are for 90 per cent confidence limit based on a A\ 2 = 2.7 criterion (Lampton, 
Margon & Bowyer 1976). 



Data (z) a E c R iV H ,G N H £ E Fc a Fc xVdof 

GO 0.017 0.90+° °l 730+^° 0.76+°^ 0.5 32±^ 530±|Jg 6.32+°^ 0.33t ?J 74.1/72 
EO 0.020 0.90±g : g? 510±||g° 0.76' 0.7 140±?tg 6.0±g-'l 0.3±g ; f 73.3/67 
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Figure 2. Comparison of the estimates of the average spectrum of RQ Scyfcrt Is obtained with various instruments. The solid-line data 
are from EXOS AT and OSSE. The solid and dashed curves are the best- fit models for the EXOSAT/OSSE and Ginga/OSSE samples, 
respectively. Note that the two curves have the actual normalizations from the fits. The agreement between the two spectra confirms 
that the X-ray variability of the AGNs is indeed fully compensated by the sizes of the samples. The dotted symbols represent the average 
HEAO-1 Al and A4 spectrum, which is in agreement with the Ginga/ EXOSAT/OSSE spectrum. The thick-line upper limits are for 
combined emission of all Seyfert galaxies observed by COMPTEL (Maisack ct al. 1995). 



spectrum is almost the same as that obtained for a larger 
sample of Seyfert Is by Maisack, Wood & Gruber (1994). 

Fig. 2 also shows the upper limits on emission of all 
Seyferts observed by COMPTEL (Maisack et al. 1995). The 
corresponding upper limits for the EXOSAT/OSSE sample 
of RQ Seyfert Is are about 3-4 times higher. We see that 
even the limits for all Sefyerts do not provide additional 
constraints on our average spectrum. The upper limits on 
7-ray emission of Seyferts were also obtained by EGRET 
for energies above 100 MeV (Lin et al. 1993). The weighted 
2-a upper limit for the EXOSAT/OSSE sample corresponds 
to EF E ^ 0.004 keV cm" 2 s" 1 at 100 MeV. This is much 
above our extrapolated models. However, that upper limit 
does rule out an extrapolation of the a — 0.9 power law 
without a break up to 100 MeV, at which energy the ex- 
trapolated power law spectrum would be a factor of ~ 5 
above the upper limit. Thus, the EGRET results provide a 
confirmation of the existence of a high-energy spectral break 
in Seyferts, independent of the OSSE results. 



5 PHYSICAL PROCESSES IN SEYFERTS 



5.1 Thermal models 

Exponentially cut-off power-laws (used here to fit the av- 
erage Seyfert- 1 spectrum) can be used to model spectra 
due to Comptonization in thermal, optically-thin, mildly- 
relativistic plasmas (e.g., Z94). The parameters, a and E c , 
can be related to the Thomson optical depth, r, and tem- 
perature, T. For E c — 400 keV, which fits both our aver- 
age spectra and IC 4329A, Z94 find E c ~ 1.6&T (implying 
kT ~ 260 keV). Z94 also provide an expression for a(r, kT) 
in a slab geometry, which yields r ~ 0.1 for a — 0.95, 
E c = 400. Thus, optically-thin, mildly relativistic plasmas 
can explain the observed spectra of Seyfert Is. The large val- 
ues of E c obtained here rule out models with Comptoniza- 
tion in optically thick plasmas (Sunyaev & Titarchuk 1980), 
used in the past to fit Seyfert-1 spectra (e.g., Miyoshi et al. 
1988). 

Haardt & Maraschi (1993) have proposed that the hot 
plasma in Seyferts forms a corona above the surface of an ac- 
cretion disk, and that most of the energy dissipation occurs 
in the corona. Then the plasma temperature can be deter- 
mined from the disk-corona energy balance, which makes 
the model more self-consistent. The disk-corona model also 
accounts for the Compton-reflection spectral components in 
the spectra of RQ Seyfert Is. Z94 found that the best-fit 
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parameters of the X7 spectrum of IC 4329A satisfy that en- 
ergy balance, i.e., the hard corona emission reprocessed by 
the disk self-consistently provides the seed of soft photons 
(in the UV range) for Compton upscattering into the hard 
spectrum. 

Those results have recently been confirmed by a more 
sophisticated treatment of the radiative transfer in Stern et 
al. (1995; hereafter S95). We apply their results on homoge- 
neous slab coronae to our average spectrum of RQ Seyfert 
Is. The spectrum has a ~ 0.9 and the overall 2-18 keV 
spectral index, CK2-18 — 0.7 (including the reflection compo- 
nent). For that ct2-i8 and assuming a pure e corona, S95 
obtain r ~ 0.05 and kT ~ 330 keV (roughly corresponding 
to E c ~ 500 keV) in agreement with E c obtained from our 
fits, see Table 1. This agreement supports the dissipative 
corona model of Haardt & Maraschi (1993). 

On the other hand, Haardt, Maraschi & Ghisellini 
(1994) point out that since the UV fluxes in many Seyfert 
Is are much larger than the X-ray fluxes (Walter & Fink 
1993) the model of Haardt & Maraschi (1993), with most 
of the dissipation occuring in the corona, is ruled out. In 
that model, the UV emission is due to reprocessing of the 
X7 emission directed towards the cold disk, and EFe in 
the UV and in X-rays are expected to be of the same order 
of magnitude. To the contrary, EF E (1375 k)/EF E {2 keV) 
[from the best-fit values in Walter & Fink (1993)] for our 
RQ sample equals 0.8, 36, 0.2, 4.2, 19, 17, and 12 for MCG 
8-11-11, NGC 3783, MCG -6-30-15, NGC 5548, ESO 141- 
G55, NGC 7469, and Mrk 509, respectively. We see that 
the ratio is > 10 for 4 objects. These large ratios can be 
explained if the corona is patchy rather than homogeneous 
and the corona dissipation dominates the disk dissipation 
only in the vicinity of an active region (a 'patch'), but not 
globally (Haardt et al. 1994). 

We also point out that a pure e pair corona would not 
form a thin slab above the disk surface (which was assumed 
by Haardt & Maraschi 1993). Hydrostatic equilibrium im- 
plies Hc/Rd = (2r0) x / 2 for a pure pair, gas pressure- 
dominated, corona, where H c is the corona scale-height, 
is the disk radius, r = R d c 2 /(2GM), and 6 ee kT/(m e c 2 ). 
Thus, H c /R,i <; 1 for Or <; 1 and the assumption of a slab 
geometry breaks down. 

S95 also provide one more argument against homoge- 
neous slab coronae based on pair equilibrium. They calcu- 
late the maximum local compactness, l\ = LiaT/(H c m e c 3 ) 7 
at which the homogeneous corona can be in pair equilibrium. 
Here L\ is the luminosity from a characteristic local volume, 
H 3 , in the corona. For 02-18 — 0.7 of the Gmga/OSSE sam- 
ple, S95 obtain £\ Ss 0.4. The maximum value is achieved in 
a pure pair corona, and a presence of ionization electrons 
lowers the equilibrium value of £\. S95 point out that this is 
in conflict with the AGN variability data (Done & Fabian 
1989). On the other hand, S95 show that higher t\ are possi- 
ble for patchy coronae. For a?2-is — 0.7 of the Ginga/OSSF, 
sample, £\ Ss 15 for an active region in the form of a hemi- 
sphere. That compactness appears compatible with the com- 
pactnesses estimated from X-ray variability (Done & Fabian 
1989). The plasma parameters of RQ Seyfert Is implied by 
thermal Comptonization are r ~ 0.2 and kT ~ 0.7 MeV 
in the case of pure pairs (see Fig. 1 in S95), which kT is 
compatible with our fitted values of E c . 

However, we point out that the results of S95 need to 



be modified to account for the difference between the lo- 
cal compactness (used to compute pair equilibria) and the 
global compactness, £ = Lot / (Rdm e (?) (where Ra is the ra- 
dius of the region where the X7 luminosity L is produced). 
The latter (rather than the former) compactness is con- 
strained from the variability data by R^ ^ cAt (Done & 
Fabian 1989). Using the formalism of Svensson & Zdziarski 
(1994) for the scale-height of the corona, assuming O ~ 0.5, 
and using the constraint on the size of the X-ray producing 
region of r ^ 20 (Tanaka et al. 1995; Fabian et al. 1994; 
Mushotzky et al. 1995), we obtain £ > 40(1 + n + /n p )~ 1/2 ^. 
Thus, I 3> t\ is possible in coronae. This largely resolves 
the conflict between the small local compactness required 
by pair equilibrium in a homogeneous corona and the vari- 
ability data. Thus, the main argument against homogeneous 
coronae remains the large UV fluxes (Haardt et al. 1994), 
rather than pair production (S95). 

More evidence against homogeneous coronae in Seyfert 
Is is, however, provided by X-ray variability. Czerny & Lehto 
(1996) find that some variability light-curves from EXOSAT 
are truly stochastic rather than due to a deterministic chaos, 
which implies X-ray emission by multiple active centers 
rather than by a single extended source. 

We note that studies of thermal pair plasmas in pair 
equilibrium predict no distinct pair annihilation even from 
a pair-dominated plasmas (Maciolek-Niedzwiecki, Zdziarski 
& Coppi 1995) . A future detection of such a feature would 
indicate the plasma is eith er nonthermal with a hard electron 
injection (see Section 5.2) or there is a strong pair wind and 



pair annihilation takes place in a region spatially different 
from the hot corona (e.g., Maciolek-Niedzwiecki et al. 1995). 



5.2 Nonthermal models 

In nonthermal models, electrons are accelerated to or in- 
jected with a nonthermal distribution extending to rela- 
tivistic energies. The electrons Compton upscatter UV seed 
photons to the X7 energy range, and the 7-rays may pro- 
duce relativistic pairs, supplementing the primary injection 
of nonthermal electrons. We use here a numerical model of 
Lightman & Zdziarski (1987) with modifications given in 
Zdziarski, Coppi & Lamb (1990). We also take into account 
Compton reflection and absorption, which best-fit parame- 
ters for models below are almost the same as for the cut-off 
power- law model, Section ^. We fit the Ginga/OSSE data 
only since they provide more stringent constraints on the 
continuum parameters than the EXOSAT/OSSE data (see 
Section §). 

The simplest nonthermal model consists of power law 
electrons with the steady-state index of p = 2a + 1 ~ 2.8 
that singly-scatter some soft seed photons in the Thomson 
regime, yielding a power-law photon distribution with a ~ 
0.9 and a high-energy cutoff at E 3> 511 keV (Blumenthal 
& Gould 1970). If the relativistic electrons lose energy in 
the Thomson regime and thermalize within the source, the 
index, p, corresponds to injection of relativistic electrons 
with an index less by one (Blumenthal & Gould 1970), V ~ 
2a ~ 1.8. Since the model gives no high-energy cutoff in the 
hundred-keV range, it fits our data worse than the thermal 
model (represented by a cut-off power law, Section ^). 

The best fit of such a model to the Ginga/OSSE data 
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shown by the solid curve in Fig. 3a corresponds to a low non- 
thermal compactness, £ = 0.1 (implying negligible e pair 
production), the compactness in blackbody photons 20 times 
the nonthermal compactness (assuring that the X7 spectrum 
is entirely due to the first-order Compton scattering), and 
the best-fit r = 1.85. We see that the model is above the 
data at ~ 150-350 keV, which results in A\ 2 = +4 with re- 
spect to the thermal model. Increasing the nonthermal com- 
pactness steepens the X7 spectrum (Svensson 1987; Light- 
man & Zdziarski 1987) and further worsens the fit. E.g., 
I = 20, T = 1.8 yield Ax 2 = +8 with respect to the thermal 
model. Thus, we reject the (r = 1.8)-injection model. 

On the other hand, Zdziarski et al. (1990) have pro- 
posed that the X-ray spectral index of a ~ 0.9-1 of Seyfert 
Is is due to a nonthermal model with dominant pair produc- 
tion. The primary electrons are injected monoenergetically, 
which implies a = 0.5 in the absence of pair production 
(Blumenthal & Gould 1970). Then, saturated nonthermal 
pair production yields the X-ray spectral index of a ~ 1 in 
the limit of I 3> 10 (Svensson 1987; Lightman & Zdziarski 
1987), which can explain the X-ray spectra of Seyfert Is. 
However, this model also predicts a steepening of the spec- 
trum at a few tens of keV due to downscattering of hard 
X-rays and soft 7-rays by thermalized, optically thick, pairs 
(Sunyaev & Titarchuk 1980). This model was ruled out for 
IC 4329A by Z94, which conclusion we extend here to the 
average Seyfert-1 spectra. The dotted curve in Fig. 3a shows 
the best fit of this model [I ~ 30). We see that the model lies 
below the ~ 50-150 keV spectrum, which results in an un- 
acceptable A% 2 = +25 with respect to the best-fit thermal 
model. 

However, we find that some nonthermal models inter- 
mediate between the two models rejected above fit our data 
satisfactorily (as also shown by Z94 for IC 4329A). We con- 
sider here models with a power-law electron injection (be- 
tween the Lorentz factors of 1.3 and 10 3 ), pair production, 
and allowing for repeated Compton scattering. This yields 
r = 5.0^3'g with the same \ 2 as f° r the thermal model. 
The spectrum of the best-fit model (I ~ 130; the solid 
curve in Fig. 36) is due to repeated Compton upscatter- 
ing by both the nonthermal and the thermal parts of the 
electron/pair distribution (Zdziarski et al. 1990; Ghisellini, 
Haardt &: Fabian 1993). Note the similarity of the spec- 
trum to that of the thermal model (Fig. 1). This nonther- 
mal model yields no distinct pair-annihilation feature. The 
annihilation feature becomes stronger with decreasing Y, il- 
lustrated by the dotted curve in Fig. 36 for V = 2.5 (£ ~ 30; 
Ax 2 = +1.6). Note that the annihilation features in Fig. 3a, 
b are allowed the OSSE data and do not constrain the mod- 
els by themselves. Decreasing V also decreases the impor- 
tance of repeated Compton upscattering. The spectrum for 
the lowest allowed F (= 1.4) is dominated by the first-order 
nonthermal Compton scattering by both the injected non- 
thermal electrons and nonthermal pairs (which pairs steepen 
the X-ray spectrum from a — T/2 — 0.7 to the observed 
value of a — 0.9). 



6 CONCLUSIONS 

We have obtained the average spectral parameters of RQ 
Seyfert Is in X-rays and 7-rays based on samples of ob- 
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Figure 3. The average spectrum of RQ Seyfert Is from Ginga 
and OSSE modeled by nonthermal plasmas, (a) Two examples 
of rejected models. The solid curve corresponds to a model with 
power-law electrons injected with the index T = 1.85 and with 
negligible both pair production and repeated Compton scatter- 
ing. The solid curve corresponds to a model with monocncrgctic 
injection of electrons and a saturated pair cascade (Zdziarski et 
al. 1990). (b) Nonthermal models with power-law injection of elec- 
trons, providing fits as good as the best-fit thermal model. The 
solid curve represents the best-fit model with T = 5.0 and negli- 
gible pair production. The dashed curve represents a model with 
T = 2.5 and moderate pair production. See Section 5.2. 



jects observed by both Ginga and OSSE and by both EX- 
OS AT and OSSE. The estimates of the average spectrum 
from both samples are fully consistent with each other. Fur- 
thermore, they are consistent with the results of HEAO-1 
Al and A4. The average spectrum contains an underlying 
power law with a ~ 0.9. There is also a Compton-reflection 
spectral component corresponding to cold matter covering a 
solid angle of ~ 1.57T as seen from the X7 source. The power 
law continues to soft 7-rays and it breaks or it is cut off with 
an e-folding energy of 250 keV. The range of the e-folding 
energy obtained for our average spectra is fully consistent 
with that of IC 4329 A (Z94; Madejski et al. 95). The aver- 
age e-folding energy is not 40-50 keV, which was reported 
earlier based on the OSSE data alone. 

Our average spectrum agrees well in X-rays with the av- 
erage spectrum of all Seyfert Is observed by Ginga (NP94). 
The intrinsic dispersion in the spectral index and the con- 
tribution of reflection in Seyfert Is is given by NP94. On 
the other hand, the limited statistics of our OSSE spectra 
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allows us to provide only the range of the e-folding energy 
for the average spectrum with a = 0.90 ± 0.05. The value 
of the e-folding energy in individual AGNs appears to cor- 
relate positively with the X-ray spectral index (Zdziarski et 
al. 1996). 

The average spectra of RQ Seyfert Is can be modeled 
by Comptonization models with either thermal or nonther- 
mal electrons. The Comptonizing plasma in thermal models 
is optically thin and relativistic. The plasma is likely to form 
a patchy corona above the surface of an accretion disk. Some 
models with relativistic nonthermal electrons predict the 
presence of an annihilation feature around 511 keV, which 
can be tested in the future by the INTEGRAL observatory 
(Winkler 1994). 
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